Resonant harmonic generation and collective spin rotations 
in electrically driven quantum dots 
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Spin rotations induced by an ac electric field in a two-electron double quantum dot are studied 
by an exact numerical solution of the time dependent Schrodinger equation in the context of recent 
electric dipole spin resonance experiments on gated nanowires. We demonstrate that the splitting 
of the main resonance line by the spin exchange coupling is accompanied by the appearance of 
fractional resonances and that both these effects are triggered by interdot tunnel coupling. We 
find that the ac driven system generates residual but distinct harmonics of the driving frequency 
which are amplified when tuned to the main transition frequency. The mechanism is universal for 
electron systems in electrically driven potentials and works also in the absence of electron-electron 
interaction or spin-orbit coupling. 

PACS numbers: 73.21. La, 03.67.Lx, 71.70.Gm, 75.70.Tj, 81.07.Ta, 42.65.Ky 



I. INTRODUCTION 

The ideai of processing quantum information stored in 
spins of electrons confined in quantum dots has motivated 
a significant theoretical and experimental effort within 
the last decade. One of the necessary pre-requisites for 
quantum gating is coherent single spin manipulation. 
Single-spin rotations can be performed using electron 
spin resonance - Rabi oscillations in external microwave 
radiation resonant with the Zeeman splitting of energy 
levels in a magnetic field {B). Electron spin resonance 
was implemented in a quantum do1»2i using an embed- 
ded on-chip microwave source. In gated quantum dots 
the microwave field has been successfully replaced by 
ac voltagesj^"— The periodic motion of the electron in- 
duced by the ac field subjects its spin to an oscillating 
momentum-dependent spin-orbit (SO) fieldj^ii^ leading 
to the electric-dipole spin resonance (EDSR)iii The spin 
rotations are detected in two-electron double quantum 
dots systems^Ti^ using the Pauli blockade of the current 
flow which occurs when the dots become occupied by elec- 
trons with parallel spins. The rotation of the spin lifts 
the Pauli blockade when the frequency of the ac electric 
flcld is tuned to the resonant transition. 

The detailed structure of the EDSR was recently 
resolved^ in a double dot produced in a gated InSb quan- 
tum wire with strong SO interactions. The experimental 
data [Fig. 2(b) of Ref. Q include a double line corre- 
sponding to transitions from the spin-polarized triplet 
r+ ground-state to a doublet formed by: 1) singlet 5, 
and 2) unpolarized triplet Tq, as well a single line at half 
the resonant frequency. Half-resonances were previously 
observed also in InAs quantum wire dots [Fig. 2(b) of 
Ref. 01 as well as in GaAs planar quantum dotS)^ Anal- 
ysis of the dipole moment induced by an ac field in the 
SOT(?Zei-subspace of two-electron systems was given in Ref. 
[12I in terms of flopping the pseudospin mode. 



In this paper we report on the solution of the time- 
dependent Schrodinger equation for the two-electron sys- 
tem in an ac fleld induced by gates in the presence of SO 
coupling. We flnd spin transitions involving both Rabi 
oscillations for the resonant driving frequency as well as 
fractional resonances which arc consistent with the exper- 
imental datai^ We show that the mechanism responsible 
for the appearance of the fractional lines is the resonant 
amplification of the higher harmonics residually present 
within the driven system. 



Solution of the time-dependent Schrodinger equation 
is one of the methodsi^"— applied in theories of high 
harmonics generation by atoms and molecules in intense 
laser fieldsi^"— in the quest for controllable sources of 
ultra energetic photons. Noble gases or simple molecu- 
lar systems (N2, O2, CO2) generate non-resonantly high 
harmonics of the driving laser field of intensity 10^^ 
W/m^ with local field amplification by plasmonic metal 
nanostructureSfiSi or 10^'^ W/m^ in standard conditions. 
The amplitude of the ac electric field applied in EDSR 
for quantum dots (a fraction of kV/cm) corresponds to a 
laser radiation of only 10^ W/m^. Nevertheless, we find a 
distinct - although residual - appearance of a second and 
third harmonics of the driving frequency ujac in the elec- 
tron motion within the double dot. We demonstrate that 
the harmonics of the driving frequency are essentially re- 
inforced when brought to resonance with the Rabi direct 
transition frequency. We indicate that this phenomenon 
is quite general for ac driven electron systems confined 
in quantum dots, in particular that it appears also for a 
single electron and in the absence of SO coupling. As a 
result, the confined system is driven into an excited state 
by frequency uj being a fraction of the excitation energy 
AE, i.e. Hlu — AE/n, which is similar to n-photon opti- 
cal transitions)^ 



2 




Fa„(t) 



' 00 


I \ 


B 


\ 



cj = —51, e = 16.5 and take the Rashba constant a = 50 
meVnm. The ac field amplitude Fac = 0.1 kV/cm is 
assumed. All the calculations are performed within the 
finite difference scheme on the {xi,X2) space with exact 
inclusion of the electron-electron correlation. 



FIG. 1. (color online) Schematic of the considered confine- 
ment potential of a nanowire double dot system. 



III. RESULTS 



II. MODEL 

The considered two-electron system is described by the 
Hamiltonian 
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with the single-electron energy operator, 
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with magnetic field B aligned along the z-dircction. The 
momentum operator is ?iki = —iKS/i as we neglect the or- 
bital effects of the magnetic field for low values of B and 
in strong confinement in the plane perpendicular to the x- 
dircction. The V{r, t) stands for the confinement poten- 
tial taken in a separable form V{r, t) = Vx{x, t)Vy^z{y, z). 
We include Rashba SO interaction^ resulting from the 
electric field - generated by the system of gates on which 
the nanowire is deposited - which is assumed perpen- 
dicular to the wire (parallel to the z-direction), Hso = 
Oi{uxky — Uykx)- Figure [T] depicts the considered confine- 
ment potential. The structure is assumed 250 nm long, 
with 

Vx{x,t) = Vsx{x) + eFbiasX + eFacxf{x) cos{ujact). (3) 

The last term in Eq. ^ represents the ac field, which 
is applied to the left dot only (see [1]), i.e. f{x) = 1 
in the left dot and outside. Vsx is a double quantum 
well potential with a 30 nm -thick barrier of height Vb 
in the center. A constant F^ias — —0.1 kV/cm is taken 
for 8 mV source-drain bias voltage.— We assume a strong 
radial parabolic confinement in the (y, z) direction which 
freezes the lateral wave functions of both electrons into 
Gaussians * = {y/nl)~^ exp [-{y^ + z'^)/2l'^] , with I = 
30 nm. Upon integration of (1) with the lateral wave 
functions one arrives at an effective Hamiltonian^i 



2m* dx' +Mx^^t)-(^'^y^kx 
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which is used in this work. Unless stated otherwise we 
apply material parameters'^ for InSb: m* = 0.014mo, 
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FIG. 2. (color online) Energy spectrum of coupled quantum 
dots with barrier height Vb — 100 meV. (b) Map of the spin 
transitions after 30 ns: the minimal value of the x component 
of the spin obtained during the simulation in h/2 units (the 
initial value is 2). (c) Evolution of the spin x component in 
the left (blue curve), right (red curve) and both dots (black 
curve) at the resonance marked with the black arrow in (a). 

Figure HJa) shows the energy spectrum of weakly cou- 
pled quantum dots (14 = 100 meV). The four- fold degen- 
eracy of the ground state at _B = is due to weak tunnel 
coupling between the dots. The degeneracy is lifted in 
nonzero B: the triplet energy levels with spin aligned 
parallel and antiparallel to the a;-direction - T_|_ and T_ 
respectively^^ (plotted with the red and blue curves) - 
are split by the Zeeman interaction. The two other states 
- To and S (plotted as the green solid and black dashed 
curves) - with zero value of spin component in the x- 
direction remain degenerate (with energy separation be- 
low 0.1 /xeV). 

We initialize system in the state. For B = 20 mT 
and ac frequency tuned to the energy difference between 
the and S states {hujac ~ Es — 55 ^eV which corre- 
sponds to the oscillation period Tac — 75 ps) wc see [Fig. 
life)] that after about 2.7 ns the spin of the electron in 
the left dot - wiggled by the ac field - is inverted, while 
the spin in the right dot remains unaffected. 

The EDSR experiments probe the spin-rotations by 
measuring the map of current leakage through the spin 
blockade as function of the driving ac frequency ujac and 
external magnetic field magnitude B. Figure [^Jb) shows 
the minimal value of the spin obtained during a time evo- 
lution of 30 ns as function of B and the driving frequency 
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(the initial x component of the spin is 2 in fi/2 units). 
A single line corresponding to the r+ — >• (5, Tg) doublet 
transition is obtained. Outside this resonant line the ac 
field does not influence the spin. Note, that the transi- 
tion T_ is not observed since it requires rotation 
of both spins. 
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FIG. 3. (color online) Same as Fig. [2]but for stronger interdot 
coupling, i.e. Vb = 17 meV. 

For stronger interdot coupling (VJ, = 17 meV) the split- 
ting of To and S energy levels (the exchange energ}iii^) 
becomes nonzero, J = Eto — Es — 5.6 /xeV [see the 
energy-spectrum in Fig. [3ja)]. For the Wac tuned to the 
T_|_ — > 5* transition we observe that the ac field applied to 
the left dot rotates the spins in both dots [Fig. Ellc)]. This 
is due to the spin-exchange interaction which is now acti- 
vated by the interdot tunneling. The characteristic spin 
swap time r = ixhl J ~ 370 ps corresponds to the inter- 
vals between the local extrema of spins observed in Fig. 
|2t^c). In the plot there are also fast oscillations of the 
spin component visible. Their period corresponds to the 
period of the ac field, i.e. t^c — 63 ps for hiOac = 66 /.teV. 
They are due to spin prccessioii^^ induced by the spatial 
electron oscillation driven due to stronger interdot tun- 
nel coupling as compared to VJ, = 100 meV case. Near 
B=0 we observe an avoided crossing of lowest-energy lev- 
els due to the SO interaction. For B > 0.1 mT they can 
be identified by their spin-x component as S, T_|_, Tq, and 
T— , see also the end of the Section. 

The map of minimal spin states encountered during 
a 30 ns simulation is presented in Fig. [3jb). At the 
diagonal of the plot two major lines emerge. They cor- 
respond to the transitions from the state to the S 
(for TiLOac = Eg) and Tq (for fioJac — Etq) states with 
spin X components. This double line was observed in the 
experimentsii^ and attributed to different gi-factors in the 
dots with the assumption that a single spin responds to 
the ac field, and the local differences in g factor are due 
to variation of the confinement composition. In fact the 
present simulation shows that the lines are split when 
J 7^ 0, which implies the coupling between the spins in 
both the dots. Rotation of the spin in the left dot to 



which the ac field is applied results in the spin rotation 
of the other electron. When only the spin in the left dot 
is inverted, the final state corresponds to a spatial "spin 
density wave" which is not an eigenstate of the spin, but 
a superposition of S and Tq states, which can be a sta- 
tionary Hamiltonian eigenstate only provided that S and 
To are degenerate as it is the case in Fig. 2. 

The transition Tj^ — > S lifts directly the Pauli block- 
ade, while the T+ — )• Tq transition lifts the blockade only 
indirectl'^ due to the mixing of S and To states - that 
arc close in energy - by the hyperfine field. The red 
curve in Fig. ^h) obtained for hui = Et_ corresponds 
to the transition to the T_ state which requires the ro- 
tation of spins in both dots, and is therefore not visible 
for Vb = 100 meV [see Fig. ^h)]. In nonzero B the T_ 
energy level is too far on the energy scale to mix with the 
S state via the nuclear spins. For that reason the transi- 
tion to T_ state does not unblock the current flow2^ and 
therefore this line is missing in the experimental datai2iS 
of the frequencies lifting the spin blockade of the current. 
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FIG. 4. (color online) (a) Spin transitions obtained for B = 
20 mT after 3 ns (black), 10 ns (green) and 30 ns (red) for 
Vi, = 17 meV (minimal x component of the spin acquired by 
the system subjected to ac field with T+ as the initial state), 
(b) Probability of transition to the 5* (black curves), To (green 
curves) and T- (blue curves) states after 30 ns. 



Besides the direct Rabi transitions additional ones for 
lower frequencies are clearly visible in Fig. [3l^b). Let us 
focus on a cross section of the map Fig. EJb) obtained for 
B = 20 mT presented in Fig. IH^a). The transition prob- 
ability is plotted in Fig. IH^b). The three broad peaks 
(marked with Es , Etq and Et^ ) correspond to direct 
Rabi transitions. The narrow resonances observed for 
lower uJac correspond to fractions of the frequencies of the 
direct transitions. The transition probabilities depicted 
in Fig. IHb) exhibit one-half and one-third (those are not 
fully saturated in the plot resolution, i.e. hl^ojac = 50 
neV) T_|. — >■ 5 transitions (black curve) and fractional 
T+ — > T_ transitions (blue curve). The fractional tran- 
sition to the 5' state (for hujac = Es/2) is of particular 
importance as it lifts the spin blockade in the EDSR ex- 
periments and this is the fractional resonant line that is 
visible in the experimental maps of Refs<^i^. The di- 
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rect Rabi oscillations are rather slow22 but the fractional 
ones are even slower. In Fig. IH^a) the transitions after 
3 ns, 10 ns and 30 ns are plotted. For 3 ns - the direct 
transitions are already fully resolved in contrast to the 
fractional ones. At the left upper corner of the map Fig. 
[DJb) one can observe additional resonance line which is a 
fractional resonance to the fourth excited state. 



10 ^ no fractional transitions are observed. 
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FIG. 5. (color online) Fourier transform of the total momen- 
tum calculated for three different driving ac frequencies hujac- 

Let us focus on the origin of the fractional transitions. 
First we consider the ac frequency hu>ac for which no tran- 
sition occurs. The Fourier transform of the total momen- 
tum is presented in Figs. [5][b,c). We observe that when 
the electron is driven by an ac field its motion is periodic 
consisting of: i) the driving frequency fiuiac, ii) its har- 
monics [marked with the red dashed lines in Fig. [S] and 
iii) the resonant frequency corresponding to the direct 
r+ — >■ S* transition [marked with the blue dashed line at 
Fig. [5]. When the driving frequency is such that one of 
its harmonics matches the resonant one [Fig. [5Ja)] its 
amplitude is greatly amplified and the system exhibits a 
resonant transition. 

We find that there is a relation between the intensity of 
the fractional line in the transition maps and the proba- 
bility of finding both electrons in the same dot in a given 
final state. In particular for VJ, = 17 meV the proba- 
bilities for subsequent states are: S - 0.02, Tq - 0.002 
and T_ - 0.04. Both the transitions for TitUac = Es/2 
and hujac = Et_/'2- occur in t ~ 25 ns with the half- 
width of the transition peak approximately 100 neV wide 
(the transition for TiLOac = is as long as 200 ns 

with the half- width of the peak about 10 neV). The frac- 
tional transitions to Tq are missing in Fig. 2] and we 
do not observe the generation of a residual frequency for 
hu! = Eto in Fig. [5] In fact the fractional transition for 
hujac = Et„/'2 does occur but with a line width that is 
narrower than the plot resolution of 50 neV - the half- 
width of the peak is about 1 neV with the transition 
time more than t ^ 1 /is. For non-zero double occupancy 
probability the electrons are at least partially driven over 
the entire double dot area by the ac field and the genera- 
tion of higher harmonics becomes effective. Consequently 
for Vb — 100 meV where the probabilities are as small as 
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FIG. 6. (color online) Transition probability to the first- 
excited (green curve) and second-excited (red curve) state 
after 10 ns for one-electron single dot (VJ, — 0) with (solid 
curves) SO coupling and without (dashed curve). Inset 
presents the single-electron energy spectrum with SO cou- 
pling included along with the direct Rabi transitions. 

The source of the fractional resonance observed in the 
experiments is the dynamics of a non-adiabatically driven 
electron system. In order to demonstrate that let us re- 
duce the problem to a single-electron one (we also lift 
the interdot barrier). The energy spectrum for such a 
system is presented in the inset to Fig. [S] The transition 
between the ground state and the first-excited state is 
only possible through a spin rotation. The transition for 
hujac ~ E2 occurs between the states of the same spin. 
In Fig. [B]we show the transition probabilities to the first- 
and second- excited states. We observe both the direct 
transitions and the fractional ones. When we switch off 
the SO coupling the transition to first excited state is 
blocked as spin becomes decoupled from the electron mo- 
tion. However the transitions to the second excited state 
along with its fractional components is still present - see 
the red dashed curve at Fig. [B] 
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FIG. 7. (color online) Same as Fig. [3] but for weaker SO 
coupling strength. The map in (b) is obtained after 60ns. 

For completeness we present the case of weaker SO 
coupling - namely we apply a ~ 25 meVnm. The energy 
spectrum for strongly coupled dots (we chose Vb = 27 
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meV to obtain similar coupling strength between S and 
Tq as previously, i.e. the exchange energy J = 6.3 ^eV) 
is presented in Fig. ^a). For lower value of a one ob- 
tains a singlet and a triplet energy levels at i? = and 
the avoided crossing between T-t- and S energy level oc- 
curs at B > 0. Weaker SO interaction results also in 
a longer spin rotation time [sec Fig. [71[c)] and now the 
T-|_ — >■ S transition occurs in 7 ns (compared to 2.7 ns 
in Fig. 3). The exchange driven small-amplitude spin 
oscillation with r = 328 ps are visible in Fig. 7(c) along 
with fast oscillation due to spin precession. In the map 
of spin transitions Fig. ^h) both the double central line 
and the half-frequency transitions to the S and r_ states 
appear. Only now the lines are narrower as compared 
to stronger SO interaction case (as the amplitude of the 
SO effective magnetic field that drives the spin transi- 
tion is decreased). The half- frequency transitions occurs 
in 54 ns (compared to 25 ns for a = 50 meV nm). Oth- 
erwise the EDSR transition map remains qualitatively 
unchanged. 

IV. CONCLUSIONS 

In conclusion we studied the electrically induced tran- 
sitions between the electron states in quantum dots. We 
presented that the electron oscillations induced by an 
ac field is accompanied by residual harmonics of the 



driving frequency. We demonstrated that the fractional 
transitions observed in EDSR experiments involve reso- 
nant amplification of the harmonics in the electron dy- 
namics when they match the Rabi transition frequency. 
Moreover we indicated that the resonant amplification of 
higher harmonics is an intrinsic phenomenon of a driven 
electron system which occurs also for a single charge and 
without SO coupling. 

In the two-electron system of the double dot when the 
ac field is only applied to one of the dots a non-zero 
interdot tunnel coupling is necessary for the fractional 
transitions to appear as it triggers the motion of both 
electrons. A consequence of the non-zero exchange en- 
ergy is the splitting of the main resonance line to Tq and 
S final states. Thus the appearance of the double res- 
onant line and the fractional resonance have a common 
origin. 
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